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ABSTRACT 

Context. The massive post-Main Sequence star W243 in the galactic starburst cluster Westerlund 1 has undergone a 
spectral transformation from a B2Ia supergiant devoid of emission features in 1981 to an A-type supergiant with a rich 
emission-line spectrum by 2002/03. The star was proposed as a Luminous Blue Variable undergoing an eruptive event 
by Clark & Negueruela (2004). 

Aims. We examine the continued evolution of W243 from 2002 to 2009 to understand its evolutionary state, current 
physical properties and the origin of its peculiar emission line spectrum. 

Methods. We used VLT/UVES and VLT/FLAMES to obtain high resolution, high signal-to-noise spectra on six epochs 
in 2003/04 (UVES) and ten epochs in 2008/09 (FLAMES). These spectra are used alongside other lower-resolution 
VLT/FLAMES and NTT/EMMI spectra to follow the evolution of W243 from 2002 to 2009. Non-LTE models are used 
to determine the physical properties of W243. 

Results. W243 displays a complex, time- varying spectrum with emission lines of Hydrogen, Helium and Lyman-a pumped 
metals, forbidden lines of Nitrogen and Iron, and a large number of absorption lines from neutral and singly-ionized 
metals. Many lines are complex emission/absorption blends, with significant spectral evolution occurring on timescales 
of just a few days. The LBV has a temperature of ~8500K (spectral type A3Ia + ), and displays signs of photospheric 
pulsations and weak episodic mass loss. Nitrogen is highly overabundant, with Carbon and Oxygen depleted, indicative 
of surface CNO-processed material and considerable previous mass-loss, although current time-averaged mass-loss rates 
are low. The emission-line spectrum forms at large radii, when material lost by the LBV in a previous mass-loss event is 
ionized by an unseen hot companion. Monitoring of the near-infrared spectrum suggests that the star has not changed 
significantly since it finished evolving to the cool state, close to the Humphreys-Davidson limit, in early 2003. 

Key words, stars: individual: W243 - stars: evolution - supergiants - stars: variables: general - winds, outflows 
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A recent addition to the catalogue of galactic LBVs is 
W243 1 (Clark & Negueruela 2004, hereafter Paper I) in the 
starburst cluster Wester lund 1 (hereafter Wd 1; Wester lund 
1961; Clark et al. 2005). Early observations included 
star G (= W243) amongst a group of OB supergiants, with 
Borgman et al. (1970) reporting to^— 15.1±0.2 2 and spec- 
tral type B0.5Ia, while Lockwood (1974) list m„~15.6; see 
alsoKoorneef (1977). Wester lund (1987) report m„~14.38 
from photometry obtained in 1960-66, but used a spectrum 
taken in 1981 to classify W243 as a B2Ia supergiant with a 
spectrum devoid of emission lines. However, Paper I found 
that by 2002/03 W243 was displaying the spectrum of an 
A2Ia supergiant, implying an apparent decrease in temper- 
ature of — 10 4 K during the intervening period, accompanied 
by the development of a rich emission line spectrum show- 
ing lines of H, He I, and Ca II. Excellent agreement was 
found between the spectral types of six other B- and A- 
type supergiants observed by both Wester lund (1987) and 
Clark et al. (2005), and Paper I concluded that W243 is 
an LBV that was quiescent when observed by Westerlund 
(1987) but has since undergone a significant outburst, mov- 
ing it redwards on the H-R diagram and (assuming constant 
bolomctric luminosity) close to the Humphreys-Davidson 
limit. 

More recently, W243 has been listed as an aperiodic 
variable with to„ — 15.73 by Bonanos (2007), while Groh 
et al. (2006, 2007) report infrared spectra with Pa7 emis- 
sion that resembles the cool phase of the LBV HR Carinae 
(Machado et al. 2002) and a K-band spectral morphol- 
ogy that is very similar to other LBVs such as AG Carinae 
(Groh et al. 2009), LBV 1806-20 (Eikenberry et al. 2004) 
and the Pistol Star (Figer et al. 1998). However, the 
strong He I emission seen at 6678A, 7065A (Paper I) and 
10830A (Groh et al. 2007) is anomalous, and Groh et 
al. (2006) note that the K-band spectrum also implies a 
higher temperature than that of a typical yellow hypergiant 
(YHG) and suggest that W243 may be evolving back to- 
wards a hotter state. Clark et al. (2008) report W243 as 
a weak (L x < 10 32 erg s _1 ) X-ray source 3 , while Dougherty 
et al. (2009) find a current mass-loss rate for W243 of 
4 — 6 x 10~ 6 M Q yr _1 from radio observations. No evidence 
is seen of the nebulosity typically surrounding LBVs (e.g. 
Weis 2003), although it is likely that such a nebula would 
be rapidly disrupted in the environment of Wd 1 . Finally, 
Ritchie et al. (2009) report evidence for photosphcric pul- 
sations in W243; this is examined further in Section 3.2.2. 

This paper presents further analysis of the Paper I 
dataset, along with analysis of five subsequent high- 
resolution Echelle spectra taken during 2004 and 
intermediate- and high- resolution spectra from 2005, 2008 
and 2009 to examine the physical state and evolution of 
W243 over a seven-year timescale. Details of the observa- 
tions and data reduction are presented in Section 2, results 
are presented in Section 3 and are discussed in Section 4. 



1 RA=16 47 07.5 <5=-45 52 28.5, J2000. 

2 We note that photometry from different studies of Wdl 
tends to show systematic differences on the order of a few tenths 
of a magnitude (Clark et al. 2009, in prep.) and the values for 
W243 listed here may not be directly comparable. 

3 This is consistent with X-ray observations of other LBVs 
(e.g. Muno et al. 2006), although some LBVs (e.g. r\ Carinae, 
HD 5980) have significantly higher X-ray luminosities (Muno et 
al. 2006; Clark et al. 2008) 



Table 1. Dates of observations, instruments and configurations 
used. 
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24/07/2008 
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14/08/2008 
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FLAMES 
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04/09/2008 


54713.0 


FLAMES 
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15/09/2008 


54724.1 


FLAMES 
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19/09/2008 


54728.1 


FLAMES 


MEDUSA, HR21 


25/09/2008 


54734.1 
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14/05/2009 


54965.2 


FLAMES 


MEDUSA, HR21 


18/05/2009 


54969.3 


FLAMES 
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a The VLT/UVES configurations list the instrument mode (RED 
or DICHR^l), central wavelength (nm) and cross-disperser 
(CD#3 or CD#4). The NTT/EMMI configurations list the ob- 
serving mode (REMD or RILD), binning and grating, while the 
VLT/FLAMES configurations list the GIRAFFE spectrograph 
mode and setup name. 



2. Observations & data reduction 

Observations of W243 were made on six nights in 2003 and 
2004 using the Ultraviolet and Visual Echelle Spectrograph 
(UVES; Dekker et al. 2000) located at the Nasmyth B 
focus of the 8.2m VLT UT2 Kueyen at Cerro Paranal, 
Chile. UVES was used in RED mode with cross-dispersers 
CD#3 (4760 - 6840A) and CD#4 (6600 - 10600 A), and in 
DICHR#1 mode with cross-disperser CD#3 (4600-6680A) 
for the red arm. All modes give R ~ 40000. Due to 
the strength of the Ha emission from W243, two R-band 
spectra were taken; a 500s integration to capture the Ha 
line, and a longer 1800s integration (2750s on 21/09/2003, 
MJD=52903.0) in which the Ha line is saturated but the 
signal-to-noise (S/N) ratio of the remaining features is im- 
proved. The S/N ratio is in excess of 100 at A > 6200 A 
for the long RED and DICHR#1 mode integrations, with 
a minimum on the UVES REDU CCD of -77 at -5850A, 
dropping below 50 at - 5500 A on the REDL CCD. The 
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Fig. 1. VLT/UVES spectrum of W243 covering 6100 - 6700A from 21/09/2003 (MJD=52903.0) with principal features discussed 
in the text marked. Rest wavelengths are indicated. Identified diffuse interstellar bands are also marked. 



S/N ratio averages ~160 for the integrations centred at 
860nm. Due to the high degree of reddening towards Wd 1 
(A v ~ 12; Clark et al. 2005) the S/N ratio of data from 
the blue arm in DICHR#1 mode (centred at 390nm) is very 
poor. The data were bias subtracted, flat fielded and wave- 
length calibrated using the UVES data reduction pipeline 4 , 
version 3.9.0, with version 4.1.0 of the Common Pipeline 
Library (CPL). 

In addition to the UVES spectra, observations were 
made on three nights in 2005, eight nights in 2008 and 
two nights in 2009 using the Fibre Large Array Multi 
Element Spectrograph (FLAMES; Pasquini et al. 2002), 
located at the Nasmyth A focus of VLT UT2. FLAMES 
was used with the GIRAFFE spectrograph in MEDUSA 
mode, with the 2005 observations using LR6 and LR8 to 
cover 6438-7184A and 8206-9400A with R ~ 8600 and 
~ 6500 respectively, and the 2008 observations using HR21 
to cover the 8484-900lA range with R ~ 16200. 600s in- 
tegrations were used, giving a S/N ratio in excess of 100. 
The FLAMES data were reduced in a similar manner to the 
UVES data, using version 2.5.3 of the FLAMES-GIRAFFE 
pipeline with individual spectra extracted using the IRAF 5 



task onedspec. Finally, we make use of ESO Multi- Mode 
Instrument (EMMI; Dckker et al. 1986) spectra cover- 
ing the 8225-8900A region from 06/2002 and 06/2003 (see 
Paper I), and a NTT/EMMI spectrum from 02/2006 cover- 
ing the 6310-7835A region. These were obtained using the 
NTT at La Silla, Chile, and data reduction is described in 
Paper I. 

A summary of all the observations is given in Table 1. 
Analysis made use of packages within the Starlink soft- 
ware suite, including FIGARO (Shortridge et al. 1997) 
and DIPSO (Howarth et al. 2003). Equivalent widths were 
measured using the DIPSO EW command and integrating 
the flux relative to a linear continuum based on estimated 
points on either side of the line. Radial velocities were calcu- 
lated by correcting for heliocentric velocity and then fitting 
Gaussian profiles to the absorption features using IRAF 
task ngaussfit. In general this produced a satisfactory fit, 
but a number of strong lines display asymmetric wings and 
the measurement of the line centre was refined by carrying 
out a second fit to the core of the line profile, excluding the 
asymmetric component of the absorption line. 



4 http:/ /www.eso.org/sci/data-processing/software/pipelines/ 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 



Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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3. Description of the spectra 

The VLT/UVES and VLT/FLAMES spectra of W243 re- 
veal a complex, time- varying spectrum containing absorp- 
tion and emission features, as well as absorption/emission 
blends with P Cygni, inverse P Cygni and double-peaked 
profiles: the region around the Ha line is plotted in Figure 1, 
with principal features marked. Over the entire spectrum, 
strong emission is seen from the Ha and H/3 Balmer series 
lines, the P&5 A10049 and Pae A9545 Paschen series lines, 
He I AA5015, 5876, 6678, 7065 and Ca II AA8498, 8542, 
8662. Weak emission is also seen from He I AA7281, 8583. 
Forbidden line emission is seen from [N II] AA5754, 6584 and 
[Fe II] A7155, while weaker [N II] A6548 and [Fe II] AA7388, 
7453 lines are also seen. A possible [Fe II] A4799 line is 
seen at the extreme blue end of the UVES spectra. Very 
weak [S II] AA6716, 6731 emission is tentatively identified. 
Many Fe II lines are seen in both emission and absorp- 
tion, with some lines displaying strong inverse P Cygni pro- 
files in the 21/09/2003 (MJD=52903.0) spectrum, a feature 
also seen in the O I A8446 line and higher Paschen series 
lines. Strong absorption lines are seen from Si II AA6347, 
6371, the O I A7774 triplet and N I multiplet 1 and 8 lines 
in the near infra-red, while many weaker absorption lines 
from singly-ionized iron group elements such as Cr II, Sc II 
and Ti II are also seen: a full list is given in Table A.l. 
Finally, the spectrum contains many interstellar absorption 
features, including diffuse interstellar bands, a very strong 
Na I AA5890, 5896 doublet and K I A7699. 

Due to the high degree of reddening towards Wd 1 
our analysis focuses predominantly on features redwards of 
~6000A where the S/N ratio is high and spectral features 
are clearly defined. The short-integration UVES RED arm, 
cross-disperser #3 spectra were used to examine the Ha 
line, which is saturated in the longer UVES integrations. 
All other spectral features in the 4760 — 6650A region were 
examined using the long-integration UVES data. At longer 
wavelengths, UVES cross-disperser #4 and FLAMES spec- 
tra were used. Unless noted in the text, rest wavelengths 
are taken from Moore (1945) or from the NIST Atomic 
Spectra Database 6 . 



3.1. Emission lines 

3.1.1. Balmer- and Paschen-series emission 

Paper I reported strengthening Ha line emission through- 
out 2001-2003, a trend that continues into early 2004. The 
Ha line, plotted in the left panel of Figure 2, displays a 
double-peaked profile in all UVES spectra, but the central 
absorption feature becomes less pronounced and the line 
strengthens predominantly in the redwards peak; by mid- 
2004 the bluewards peak is at virtually the same level as 
the 2003 spectrum. Broad electron-scattering wings are vis- 
ible, extending to ±1000kms _1 . The H/3 line (not shown) 
also strengthens in early 2004 and displays a similar double- 
peaked profile with a central absorption feature blueshifted 
50kms _1 from its rest wavelength. Table 2 lists the equiva- 
lent width of the Ha line from each observation, along with 
the radial velocities of the two peaks, the central absorption 
feature and the separation between the two peaks in the 
profile. Lower resolution VLT/FLAMES spectra obtained 
in 2005 show a weaker but still clearly asymmetric Ha line 




52903.0 

53098.4 

53165.2 

53196.2 
53276.0 



6554 6556 6558 6560 6562 6564 6566 
wavelength 



10046 1 0048 
wavelength 



Fig. 2. Evolution of the Ha (UVES + FLAMES, left panel) 
and Pa5 A10049 (UVES, right panel) emission lines. Insets show 
magnified views of the centre of the emission line, and rest wave- 
lengths are marked in both panels. 




http:/ /physics. nist.gov/PhysRefData/ASD/lines_form. html 



Fig. 3. The Pall A8863 feature over the course of the 
VLT/FLAMES observations. The left panel plots the region 
around the line centre (with velocities relative to the Pall rest 
wavelength) with spectra offset vertically according to the date 
of the observation. The right panel overplots four spectra (three 
from 2008, one from 2009), with the rest wavelength indicated. 
The radial velocity of the Pa<5 emission line in the 21/09/2003 
(MJD=52903.2) VLT/UVES spectrum is marked with an ar- 
row, showing its close correspondence with the emission feature 
in the Pall line and lack of long-term variability. 



profile but cannot resolve the central absorption feature 
seen in the UVES spectra, while by 2006 a NTT/EMMI 
spectrum shows that the Ha line has weakened further to 
2002/03 levels; in this case the resolution is insufficient to 
see if the asymmetry is still present. 

Paschen series lines are also seen in emission, with 
a strong Y&5 A10049 line (shown in the right panel of 
Figure 2) and weakening emission from the Pa8 A9546 
and Pa9 A9229 lines. Unlike the Balmer-series lines, the 
Paschen-series lines do not show an obvious double-peaked 
profile and the centre of the Pa<5 A10049 line is almost 
constant between observations; values are again listed in 
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Table 2. Equivalent widths (A) and radial velocities (kms -1 ) for the Ha, Pa<5 and He I A7065 emission lines"' 6 . 



Date 


EW 


r V / 


net 

v r 




AV r 


I CLU 

EW C 


v r 


He I A7065 
EW V r 


21/09/2003 


-22.3±0.2 


-75.4 


-51.6 


-26.0 


49.4 


-1.10±0.08 


-52.7 


-0.79±0.03 


-48.9 


03/04/2004 


-25.4±0.2 


-71.2 


-58.4 


-32.8 


38.4 


-1.13±0.08 


-50.2 


-0.49±0.05 


-35.0 


09/06/2004 


-25.2±0.2 


-73.8 


-61.4 


-34.0 


39.8 


-1.15±0.08 


-51.7 


-0.86±0.03 


-45.6 


10/07/2004 


-24.5±0.2 


-74.2 


-59.2 


-31.3 


42.9 


-1.03±0.08 


-50.2 


-0.87±0.03 


-43.3 


28/09/2004 


-23.9±0.2 


-72.6 


-61.6 


-31.1 


41.6 


-1.05±0.08 


-49.3 


-0.84±0.03 


-44.1 


29/05/2005 


-21.2±0.3 






-34.3±6 






-0.80±0.06 


-47.1 



°For the double-peaked Ha line radial velocities are given for the bluewards and redwards peaks, V r (b) and V r (r), the central 
absorption feature, V r , and the separation of the two peaks, AV r . 

b Except where noted, errors are conservatively estimated at ±4kms _1 for the Ha radial velocities and < ±8kms _1 for the 
He I A7065 and Pa<5 radial velocities. 

c The equivalent width of the Pa<5 line includes the absorption components seen in Figure 2 on either side of the emission line. If 
these are excluded, then the equivalent width is broadly constant at ^1.35±0.15 A (with increased errors due to the difficulty of 
estimating the continuum). 



Tabic 2. The line strengthens and narrows slightly dur- 
ing 2004, but the changes in equivalent width are predom- 
inantly due to the variable absorption on both sides of 
the emission line. This is redshifted relative to the emis- 
sion line in the 2003 spectrum, but is increasingly seen on 
the blue side of the line during 2004, with a pronounced 
P Cygni profile visible in the 10/07/2004 (MJD=53196.2) 
spectrum. The higher Paschen-series lines are seen predomi- 
nantly in absorption, with the emission that dominates the 
Pa<5 line seen instead as weaker core emission that splits 
the blue wing of the absorption line; this can be seen in 
Figure 3, which plots the Pall line over the course of the 
VLT/FLAMES observations. The absorption component 
shows notable variability (discussed further in Section 3.2) 
while the emission component is largely static, showing lit- 
tle change in strength or radial velocity over the course of 
our VLT/UVES and VLT/FLAMES spectra. 



3.1.2. He I emission 

He I AA5015, 5876, 6678, 7065, 7281 and 8583 are seen 
in emission, with strong infra-red He I A1.083^m and 
He I A2.058^m emission lines also reported (Groh et al. 
2006, 2007). No He II features are detected. The triplet 
He I A7065 line is plotted in Figure 4, with measurements 
of the equivalent width and radial velocity given in Table 2. 
With the exception of the 21/09/2003 (MJD=52903.0) 
spectrum, the line shows a near-constant profile on the 
red side of the emission line with emission extending 
~100kms _1 redwards of the line centre. The bluewards 
emission is more variable, decreasing substantially in the 
03/04/2004 (MJD=53098.4) spectrum before recovering to 
earlier levels; this decrease, which clips the peak emission, 
is responsible for the decrease in measured radial velocity 
listed in Table 2. For other epochs, the measured radial ve- 
locity is very similar to that measured from the PaS line, 
with any differences likely due to the blue-wing absorption. 
The triplet He I A5876 line, shown in Figure 5, displays 
significant variability, developing a strong P Cygni profile 
in the 03/04/2004 (MJD=53098.4) spectrum and a double- 
peaked profile separated by ^80kms _1 in the 09/06/2004 
(MJD=53165.2) and 10/07/2004 (TJD=53196.2) spectra. 
The singlet He I AA5015, 6678 lines are also variable, al- 
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Fig. 4. Evolution of the He I A7065 line over the course of the 
UVES observations. Spectra are offset vertically in the left panel, 
and overplotted in the right panel. Velocities are relative to the 
rest wavelengths of the lines, and the radial velocity of the PaS 
line in the 21/09/2003 (MJD=52903.2) UVES spectrum is indi- 
cated with an arrow. 



though the double-peaked profile is not seen and the lines 
instead weaken to barely detectable levels before recovering 
in the final spectrum. 

The lowest He I triplet term, 2s 3 S, is metastable, and 
n(2s 3 S) becomes large as the level is populated by the 
recombination cascade from Hc + (Ostcrbrock & Ferland 
2006). From here, collisional excitation populates the 2p 3 P° 
upper energy level of the He I A1.083^m line as well 
as the higher triplet levels (including the upper levels of 
He I AA5876, 7065) and the singlet levels responsible for 
the observed He I AA5016, 6678, 2.058/im emission ( Bray 
et al. 2000; Osterbrock & Ferland 2006) 7 . The 2p 3 P° level 
is the lower level of both the He I A5876 and A7065 lines, but 
the oscillator strength of the He I A5876 line is almost an 



7 Absorption of A3889 photons by the 2s 3 S level also popu- 
lates the 3p 3 P° level, with ~10% of decays proceeding via the 
3s 3 S level to 2p 3 P°, strengthening the He I A7065 emission line, 
van Helden (1972) also proposes that the 3p P° level may be 
populated via O II A539.1 fluorescence from the He I Is S level, 
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Fig. 5. Evolution of the He I A5876 line over the course of the 
UVES observations. Spectra are shown in the same manner as 
Figure 4. 
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Fig. 6. Comparison of the He I AA5876, 6678, 7065 emis- 
sion lines on 03/04/04 (TJD=13099; left panel) and 09/06/04 
(TJD=13166; right panel). The radial velocity of the Pa<5 line 
in the 21/09/2003 (MJD=52903.2) UVES spectrum is indicated 
with an arrow. 



order of magnitude greater and absorption is correspond- 
ingly more pronounced; this is clearly visible in Figure 6, 
which shows the differing absorption bluewards of the two 
emission lines. Finally we note that the weak He I A8583 
emission arises from a transition from the 10d 3 D level at 
-24.45cV to the 3p 3 P° level (van Helden 1972). Other He I 
emission lines from the high singlet and triplet levels would 
be expected, e.g. the related He I A8777 line, but, given their 
low intensity, detection requires a combination of low-noisc 
spectra and separation from other photosphcric, interstellar 
and telluric features, and none are observed. 
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Fig. 7. The evolution of the complex Ca II A8662/Pal3 feaure 
in the VLT/FLAMES spectra of W243. The left panel plots the 
region around the Ca II A8662 line (with velocities relative to 
the rest wavelength), offset vertically according to the date of 
the observation. The rest wavelength of the Pal3 line is also in- 
dicated, and the radial velocity of the Pa<5 line in the 21/09/2003 
(MJD=52903.2) UVES spectrum relative to the Ca II and Pal3 
rest wavelengths is indicated with arrows. The right panel over- 
plots four spectra (three from 2008, one from 2009) with rest 
wavelengths indicated. 



3.1.3. Fe II, Ca II and Mg I emission 

Figure 7 shows the complex Ca II A8662/Pal3 blends over 
the course of ten VLT/FLAMES spectra of W243. The 
Ca II emission comes from the multiplct-2 4 2 P°-3 2 D tran- 
sition, the upper levels of which are fed by the Calcium-H 
and -K absorption lines (which at AA3933, 3968 are beyond 
the coverage of our spectra). Unlike the YHG IRC +10 420 
(Oudmaijer 1998; Humphreys et al. 2002) we do not ob- 
serve the accompanying [Ca II] AA7292, 7324 (multiplct IF) 
emission that would result from a subsequent transition to 
the ground state; this is possibly masked by heavy telluric 
contamination, but is more probably a result of collisional 
de-excitation. In the majority of spectra the peak Ca II 
emission again falls at the same radial velocity as the Pa5 
and He I emission, blueshifted with respect to a Ca II ab- 
sorption feature that is itself heavily blended with the emis- 
sion component of the adjacent Pal3 line. However, unlike 
the Paschen-series lines, the Ca II emission extends blue- 
wards >100kms _1 from the line centre and significant vari- 
ability can be seen over short timescales in the 2008 dataset, 
while the central Ca II absorption is entirely absent in the 
2009 spectra. Weak Mg I A8807 and very weak Mg I A8837 
emission is also seen. These lines arise from upper energy 
levels at 5.75eV and 7.52eV, near the ionization potential of 
Mg I, and are due to recombination of Mg + and the resul- 
tant cascade to the ground state. Both the Mg I and Ca II 
lines are likely to be formed in a neutral H I region; the 
ionization potentials of Ca II (11.9eV) and Mg I (7.65eV) 
are too low for significant recombination to take place in 
an H II region. 

Finally, a large number of Fc II lines arc also visible 
in Figure 1. The strongest Fe II features are lines from 



which would provide a second mechanism for strengthening the 
He I A7065 emission. 



B. W. Ritchie et al.: The LBV Wd 1-243 



7 




x 1.2 



2903.' 
53098.4 - 
53165.2 - 
53196.2 
53276.0 




52903.0 




53098.4 




53165.2 


i 


53196.2 




53276.0 












Iff 









6582 6583 6584 

wavelength 



9994 9995 9996 9997 9998 
wavelength 



Fig. 8. Evolution of the Fe II multiplet-40 and -74 lines over the Fig. 10. [N II] A6583 (left panel) and the Lya-pumped 
course of the UVES spectra. Spectra are labeled with Modified Fell A9998 (right panel) emission lines in the UVES spectra 
Julian Dates, and the rest wavelengths of the principal Fe II and of W234. Rest wavelengths are marked, and the location of the 
N I lines are shown. Pa<5 line in the 21/09/2003 (MJD=52903.2) spectrum is indi- 

cated with an arrow. 
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Fig. 9. Comparison of the Fe II A6456.4 (multiplet 74, left 
panel) and Fe II A6516.1 (multiplet 40, right panel) emis- 
sion/absorption blends. Two epochs are shown, 21/09/2003 
(MJD=52903.0, solid line) and 10/07/2004 (MJD=53196.2, 
dashed line). Rest wavelengths are marked, and the location of 
the Pa<5 line in the 21/09/2003 (MJD=52903.2) spectrum is in- 
dicated with arrows. Note the similarity between the behaviour 
of the Fe II A6546.4 line to the Pall line plotted in Figure 3. 



multiplets 40 (a 6 S-z 6 D, Supper ~4.8eV) and 74 (b 4 D-z 4 P, 
Xupper ~5.85eV), which display prominent inverse P Cygni 
profiles in the initial 21/09/2003 (MJD=52903.0) spectrum 
with emission components again displaying near-identical 
radial velocities to the other emission lines in the spec- 
trum of W243. These Fe II lines show significant varia- 
tions in strength and morphology; a number of these lines 
are plotted in Figures 8 and 9, which show their evolution 
over the course of the UVES spectra. The initial inverse 
P Cygni profiles become less pronounced and the lines of 
multiplet 40 are seen predominantly in emission, displaying 
notable P Cygni profiles in the 10/07/2004 (TJD=53196.2) 
spectrum before the inverse P Cygni profiles reform. In con- 
trast, the multiplet 74 lines arc seen mainly in absorption, 



with weak emission components on the blue wing that are 
very similar to the higher Paschen-series lines plotted in 
Figure 3. A number of weak Fe II lines display no obvious 
emission (e.g. Fe II AA6175, 6332, multiplets 200 and 199 
respectively). These lines all have relatively high upper en- 
ergy levels, e.g. x 4 F and x 4 G at ~8.2eV, multiplets 199 and 
200, or the unclassified multiplet containing Fe II A6644 
with an upper level at ~9.7eV. At shorter wavelengths there 
are emission/absorption blends due to Fe II AA4924, 5018, 
5169 (multiplet 42) and Fe II AA5235, 5276, 5316 (multiplet 
49) , while a large number of weaker Fe II absorption lines 
below 6000A from multiplets 41, 46-48 and 51-57 cannot 
be examined in detail due to the poor S/N ratio. Notably, 
no emission or absorption lines of either Fe I or Fe III are 
observed in the spectra of W243, implying that the iron 
must be almost entirely in the Fe + state. 

3.1.4. [N II] and [Fe II] emission 

The spectrum of W243 displays [N II] AA5754, 6548, 
6583 forbidden lines (the latter plotted in Figure 10), 
with the [N II] A6548 line a factor of ^3 weaker than 
[N II] A6583 as expected from the relative radiative tran- 
sition probabilities. The [N II] lines display virtually un- 
changing line profiles with a FWHM ^25kms _1 and near- 
static line centres that are blueshifted by ~50kms _1 from 
rest wavelengths. [Fe II] A7155 is also seen, while weaker 
[Fe II] AA7388, 7453 lines are observed. At the extreme blue 
end of the 21/09/2003 (MJD=52903.0) spectrum possible 
[Fe II] A4799 emission is detected, but this is not observed 
in subsequent (lower S/N ratio) spectra. In the first and 
last UVES spectra the [Fe II] A7155 line appears very sim- 
ilar to the [N II] emission with virtually identical FWHM 
and radial velocities, but the intermediate spectra display 
more complex, distorted profiles. Very weak [S II] AA6716, 
6731 lines are tentatively identified. The [Fe II] A7155, 7388, 
7453 lines come from the multiplet 14F a 4 F-a 2 G transition 
(Moore 1945) with an upper energy level of 1.96eV, and 
the tentatively-identified [Fe II] A4799 line arises from the 
multiplet 4F a 6 D-b 4 P transition from a level at ~2.7eV to 
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the ground state. No higher-excitation lines (e.g. [Fe III], 
[S III]) are observed, implying a temperature of ~10 4 K in 
the line-forming region (e.g. Lamers 1988). 

3.1.5. Lya and Ly/3 fluorescence emission 

In addition to the Fe II lines from multiplcts 40 and 74, 
notable near infra-red emission lines are seen in the 
21/09/2003 (MJD=52903.0) spectrum from Fc II AA8451, 
8490, 8927, 9998 and O I A8446. These arise from fluores- 
cence by Lya and Ly/3 photons (Damineli 2001), which 
(like the He I emission lines) implies the presence of a 
hot, ionizing source. The Fe II A9998 line (plotted in the 
right panel of Figure 10) is strong in all UVES spectra, dis- 
playing moderate variability including a notable strength- 
ening accompanied by a shift redwards in the 03/04/2004 
(MJD=53098.4) spectrum. This line arises from Lya pump- 
ing from the a 4 G level (notably the lower level of the 
strong multiplet 46 and 49 lines seen in the spectrum of 
W243) to an upper energy level at ~13.4eV, which subse- 
quently feeds the b 4 G upper level of the Fe II A9998 line 
(Johansson 1984). A weaker Fe II A9956 line, produced 
by the same mechanism, is tentatively identified. Sigut & 
Pradhan (1998) predict other fluorescence lines in the 8500- 
9200A region, resulting from Lya fluorescence from the a 4 D 
level. These include Fe II AA8490, 8927 which are tran- 
sitions from Lya-pumped 4 F and 6 F levels at ~11.4eV: 
both lines are observed in the spectrum of W243. Other 
predicted lines in the 9100-9200A range fall in a region of 
our spectra that is heavily contaminated by telluric fea- 
tures and are not observed. We also observe weak emis- 
sion from Fe II A8451, which shares a common upper level 
with Fc II A8490 (Damineli 2001), although this is blended 
with a telluric absorption line. The majority of the fluo- 
rescence lines are only clearly detected in the 21/09/2003 
(MJD=52903.0) spectrum, with only Fc II A9998 seen 
strongly in all UVES spectra and Fe II A8927 seen in the 
UVES and 2008 FLAMES spectra but absent in the 2009 
FLAMES spectra. The other fluorescence lines are either 
undetectable or, at best, tentatively detectable in the in- 
termediate UVES spectra, and are only weakly observed in 
the final spectrum. 

The O I A8446 line displays an inverse P Cygni profile 
in the 21/09/2003 (MJD=52903.0) spectrum that is very 
similar to that seen in the higher Paschen-series and Fe II 
multiplet- 74 lines. The O I A7774 triplet is in absorption 
and the O I A9256 triplet, although heavily blended with a 
telluric feature, also appears to be in absorption. This sug- 
gests that the O I A8446 emission component is a result of 
Ly/3 fluorescence (Bowen 1947); the pumped 3d 3 D° level 
is linked to the upper level of O I A8446 via the permitted 
O I A11286 line (outside the coverage of both our spectra 
and the spectrum of Groh et al. 2007), whereas the 3d 3 D° 
level is linked with the upper level of the O I A7774 triplet 
via the [O I] A9205 line, which is not observed. Ionization 
balance calculations (Grandi 1980) show that Ly/3 fluores- 
cence of O I cannot occur in either an H I region (due to 
the lack of Ly/3 photons) or an H II region (as the fluo- 
rescence mechanism requires oxygen to be predominantly 
neutral, and H I and O I have very similar ionization po- 
tentials). Instead, the Ly/3 pumping of O I was found to 
only be effective in dense, partially-ionized transition zones 
between H I and H II regions (Kwan 1984). The O I A8446 



line is not seen in emission or absorption in the spectrum 
of 03/04/2004 (MJD=53098.4), recovering predominantly 
in absorption in subsequent spectra in a similar manner to 
the Fe II multiplet 74 lines discussed in Section 3.1.3. 

Finally, we note that Groh et al. (2006) report 
Mg II A2.14/mi emission in their spectrum of W243. This 
is another Ly/3 fluorescence line, with pumping from the 
Mg + ground state populating the 5p 2 P° upper level of the 
observed Mg II line (Hamann & Simon 1988). 

3.2. Absorption lines 

3.2.1. Neutral and singly-ionized metals 

VLT/FLAMES spectra of W243 from 2008 and 2009 cov- 
ering 8500-8950A are plotted in Figure 11, along with 
a VLT/UVES spectrum from 2004 covering 8670-8950A. 
Strong N I absorption lines arc seen from the 3s 4 P-3p 4 D° 
(multiplet 1) and 3s 2 P-3p 2 P° (multiplet 8) transitions, 
along with strong Paschen-series Pall . . . Pal6 absorp- 
tion lines. Ca II AA8498, 8542, 8662, 8927 display emis- 
sion/absorption blends, although the Ca II A8912 line is 
seen in absorption only; this line originates from the same 
multiplet as Ca II A8927, and it is likely that the emission 
feature in the latter is a blend with a Fe II A8927 emission 
line (e 4 D-5p 4 D) that forms part of a recombination cascade 
from a Lya-pumped level at ^10eV. Notably, this emission 
component is absent from the Ca II A8927 line in the 2009 
FLAMES spectra. 

The strong near-IR N I absorption lines and lack of N II 
features in the VLT/UVES spectra of W243 imply a spec- 
tral type no earlier than AOIa (Munari & Tomasella 1999, 
see also the disappearance of N II in spectra of the LBV 
HR Car as it cools from 20kK to lOkK; Nota et al. 1997 and 
Machado et al. 2002). This is consistent with the strong 
Si II AA6347, 6371 absorption lines visible in Figure 1, which 
are are a good indicator of temperature for cool stars. These 
lines are strongest at around lOkK (i.e. spectral type ~A0), 
weakening as the temperature rises (Nota et al. 1997; 
Davies et al. 2005). Ne I, an additional indicator of tem- 
pratures greater than ^10kK, is also absent. Comparison of 
the N I line strength in the VLT/UVES and VLT/FLAMES 
spectra in Figure 11 (along with VLT/FLAMES spectra ob- 
tained in 2005 and NTT/EMMI spectra from 2003, neither 
shown) and Si II line strengths from the VLT/UVES and 
2006 NTT/EMMI spectra suggest that W243 has shown 
little change in state since it finished evolving to the cool 
phase in 2002-3 (Paper I), although we cannot rule out more 
significant variability during the periods where we do not 
have spectroscopic coverage. 

3.2.2. Pulsations 

Figure 12 shows the N I A8711 absorption line from the 
2008 and 2009 VLT/FLAMES datasets, with spectra off- 
set vertically according to their acquisition date to high- 
light changes in the line profiles. The line displays signif- 
icant variability in strength, profile and line centre, and 
similar variations are seen in the Pall line, plotted in 
Figure 3, the Pal3 line visible in Figure 7 and the strong 
Si II AA6347,6371 doublet and O I A7774 triplet lines (not 
shown). Over the course of sixteen epochs of VLT/UVES 
and VLT/FLAMES data the measured radial velocities 
from the strong neutral/singly-ionizcd metal lines lie within 
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Fig. 11. VLT/FLAMES (top: 18/5/2009, MJD=54969.3, middle: 29/6/2008, MJD=54646.2) and VLT/UVES (bottom: 28/9/2004, 
MJD=53276.0) spectra of W243. The rest wavelengths of principal features discussed in the text are indicated, and identified diffuse 
interstellar bands are also marked. 



a range of -43kms 



-18kms 1 (errors < ±3kms 1 ) with 



individual lines consistent to within a few kms in each 
spectrum; only the low-excitation Fe II multiplcts display- 
ing P Cygni and inverse P Cygni profiles in Figure 8 are 
discrepant, displaying red-shifted absorption components 
relative to the other metal lines, most probably as a result 
of significant infilling from the blue- wing emission. 

As noted in Ritchie et al. (2009), the short-term radial 
velocity variations are hard to reconcile with orbital mo- 
tion, as the implied orbital period is at odds with the nar- 
row range of measured radial velocities unless we are view- 
ing the system at sin i^0. In addition, the clear changes 
in line profile and strength apparent in Figure 12 are 
hard to understand in this scenario, and we instead inter- 
pret the radial velocity variations as resulting from bulk 
motions in a pulsating, dynamically-unstable photosphere. 
Such pulsations are well-documented in YHGs, for example 
p Cassiopaeia (Lobel et al. 1998), and very similar changes 
in N I and Fe I absorption lines are seen in VLT/FLAMES 
spectra of the pulsating YHG W265 (Ritchie et al. 2009). 
The limited spectral coverage and lack of contemporane- 
ous photometry makes determination of a pulsation pe- 
riod for W243 impossible; the behaviour of the N I lines 
in the VLT/FLAMES spectra is consistent with a long 
pulsational period similar to that observed in YHGs (de 



Jagcr 1998), but the photospheric lines also display rapid 
changes in both radial velocity and profile over a matter of 
days (e.g. the change between 18/07/2008, MJD=54665.0 
and 24/07/2008, MJD=54671.1 visible in Figures 3 and 12) 
that could imply that we are observing the effects of sparse 
sampling of a more rapid pulsational period. 

A further similarity with the YHGs p Cas and W265 
comes from the weakening and broadening of the metal ab- 
sorption lines through the pulsation period, and in particu- 
lar the development of the broad blue wing seen in the near- 
IR N I lines on 10/07/2004 (MJD=53196.2, see Figure 13); 
the O I and Si II absorption lines display a similar effect. 
The observed line broadening appears periodic, and is also 
seen in the N I lines in the 18/05/2009 (MJD=54969.3) 
VLT/FLAMES spectrum plotted in Figure 11, most clearly 
in the N I and Paschen-series absorption lines and the sig- 
nificant reduction in core emission in the N I triplet at 
^8680A. In p Cas this behaviour is correlated with max- 
imum effective temperature (Lobel et al. 1998), and re- 
sults from absorption in the wind due to a periodic phase 
of enhanced mass- loss 8 . W243 therefore appears to be in a 
relatively quiescent pulsational state similar to warm-phase 



8 Note that this is not the strongly enhanced mass loss phase 
encountered by YHGs on timescales of a few decades, accompa- 
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Table 3. Model parameters for W243. 
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Fig. 12. Evolution of the N I A8711 absorption lines during the 
VLT/FLAMES observations. As in Figures 3 and 7, the left 
panel plots the region around the line centre, with velocities rela- 
tive to the rest wavelength, for all VLT/FLAMES spectra, offset 
vertically according to the date of observation. The variability 
in line centre and profile can be clearly seen, with a broad blue 
wing extending to >100kms _1 in the 2009 observations. The 
right panel overplots four spectra (three from 2008, one from 
2009) with the rest wavelength indicated. 



Parameter 


Value" 


L, 


7 3 x 10 5 fd/4 5kDc) 2 L^ 


R, 


376 9 fd/4 SkDclR^ 


R2/3 


450.0 (d/4.5kpc)R 


T a 
1 eff 




M 


6.1 X 1U ( i a/4.okpcJ M yr 




1 n cr 1 — 1 

lbokms 


v sini 


lOkms -1 


log g 


~0.65 c 


H/He 


5 


N/N 


12.1 


o/o 


0.11 


Ca/Ca 


2.49 


Mg/Mg 


1.34 


Si/Si© 


1.51 


Fc/Fe 


1.22 



"Values for Voo, v sini and H/He are assumed, all others are 
derived. Values are for a distance of 4.5kpc (Clark et al. 2005; 
Crowther et al. 2006), with scaling listed where appropriate. 
''Limits are 8300< T c s <9300K, favouring a value at the lower 
end of this range. 

c Best fit to the higher Paschen-series absorption lines. 




Fig. 13. UVES observations of the N I absorption line com- 
plex around ~8700A. on 10/07/2004 (MJD=53196.2, solid line) 
with preceeding (09/06/2004, MJD=53165.2) and subseqent 
(28/09/2004, MJD=53276.0) spectra plotted with dashed lines 
for comparison. Rest wavelengths are marked. 



YHGs, with possible periods of increased mass-loss but no 
significant changes in state apparent during our monitoring. 

3.2.3. Modeling 

Modeling of W243 was carried out using the non-LTE 
line-blanketed radiative transfer code CMFGEN (Hillier & 
Miller 1998a; Hillier & Miller 1999), which computes line 
and continuum formation under the assumptions of spher- 
ical symmetry and a steady-state outflow. Examples of the 

nied by a significant movement redwards on the HRD (de Jager 
1998;Lobel et al. 2003) 



use of CMFGEN to model LBVs and detailed discussion 
of the workings of the code can be found in Najarro et 
al. (2009) and Groh et al. (2009) and extensive refer- 
ences therein. For W243, the H/He abundance was fixed at 
5, while C was assumed to be depleted due to the lack of 
observed lines. Ti, Cr and Ni have been assumed to have 
solar abundance, whereas abundance determinations for Fe, 
Ca, Mg, Si, N and O are firmly derived. We find that the 
absorption-line spectrum of W243 is well fit by a cool hy- 
pergiant model with a weak, unclumped wind 9 and low 
mass-loss rate. Model parameters are given in Table 3; a 
distance of 4.5kpc (Clark et al. 2005; Crowther et al. 2006) 
is assumed, with distance-scaling relationships given as ap- 
propriate (Hillier et al. 1998b). The value of log g is uncer- 
tain, with the listed value of ~0.65 providing the best fit to 
the relatively uncontaminated higher Paschen-series lines 
bluewards of 8500A. Fi gurc 14 shows the model fit to the 
strong N I multiplet 1 lines and the adjacent Pal2 A8750 
line (blended with N 1 A8747), while Figure 15 shows the fit 
to the strong Si H A6347, 6371 doublet and O 1 A7774 triplet 
lines. All of the strong metal absorption lines are well fit by 
the model, although the Pal 2 line shows significant infill- 
ing from the emission component discussed in Section 3.1.1. 
Figure 16 shows three strong Fe H multiplet 40 and 74 lines 
(two weaker lines from multiplet 199 and a multiplet un- 
classified by Moore (1945) also visible, along with several 
weak DIBs), along with the Ca II AA8912, 8927 doublet 
and overlapping Fe II A8927 emission line. In both cases 
the model reproduces the absorption component well, but 
does not reproduce the emission components. Similarly, the 
model predicts Paschen- and Balmer-series lines also in ab- 
sorption and very weak He I absorption features, rather 
than the observed emission. Along with the Lya-pumped 
Fe II and O I emission lines, this implies that the emission- 

9 For the velocity law we assume V 00 = 165kms~ 1 , appropriate 
for a quiescent YHG, and /3=1 
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Fig. 14. Comparison of the model fits to the strong N I mul- 
tiplet 1 absorption lines and adjacent Pall A8750 line. The 
VLT/FLAMES spectrum is from 04/09/2008 (MJD=54713.0) 
and rest wavelengths are marked. 
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Fig. 15. Model fits to the Si II AA6347, 6371 doublet (left panel) 
and O I A7774 triplet (right panel). The VLT/UVES spectrum 
is from 21/09/2003 (MJD=52903.0) and rest wavelengths are 
marked. 



line spectrum does not originate with the cool star but in- 
stead has a secondary source; this is discussed further in 
Section 4. 

A strong lower limit on temperature come from the 
Mg I A8807 line, while upper limits come from the He I and 
N I lines (the strength of the latter also used as a tempera- 
ture diagnostic in Paper I) . Precise temperature determina- 
tion is complicated by the presence of emission components 
in both He I and Mg I, and we find 8300< T eS <9300K, 
(i.e. spectral type A2-4Ia + ), favouring a value of ~8500K 
(A3Ia + ). At lower temperatures the Mg I absorption com- 
ponent would become too strong, even considering infill- 
ing from the superimposed emission visible in Figure 11, 
while at higher temperatures the unblended N I absorption 
lines become too weak. N is highly overabundant (12.1 x so- 
lar), with O depleted (0.11 x solar) and C highly depleted. 
This is clearly indicative of CNO-processed material being 
present at the surface of W243, with abundances similar to 




6430 6440 
Wavelength 



8920 8930 
Wavelength 



Fig. 16. Example of model fits to emission/absorption blends. 
The left panel shows three strong Fe II multiplet 40 and 74 
lines, while the right panel shows the Ca II AA8912, 8927 
doublet and superimposed Fe II A8927 emission line. In 
both cases the absorption component is well fit by the 
model, but the cool star cannot reproduce the emission 
component. The VLT/UVES spectrum is from 21/09/2003 
(MJD=52903.0) and the VLT/FLAMES spectrum is from 
04/09/2008 (MJD=54713.0). Rest wavelengths are marked, and 
the Fe II multiplets are also indicated in the left panel. 



those determined for AG Car (Groh et al. 2009). Ca is also 
notably overabundant, with Mg, Si and Fe all moderately 
enhanced relative to solar values. 



4. Discussion and Conclusions 

4.1. The spectrum of W243 

We summarize the spectrum of W243 as follows: 

- H, He I, O I, Mg I, Ca II, Fe II, [N II] and 

[Fe II] are seen in emission. Strong Hcv, H/3 and P&5 
emission is seen, with the Hcv and H/3 lines displaying 
double-peaked profiles, while the higher Paschen-series 
lines show core emission on the blue absorption wing. 
He I emission is highly variable, while O I A8446 and 
many Fe II lines also show significant variability. [N II] 
and a few [Fe II] lines are observed. Apart from a few 
highly-distorted lines, all permitted and forbidden emis- 
sion lines display a blueshift of ^50kms _1 with negligi- 
ble epoch-to-epoch variability. 

— Many neutral or singly-ionized metals are seen 
in absorption. Strong absorption lines of Si II, N I 
and O I are seen, along with many Fe II lines; Fe I 
and Fe III features are not observed. Many weaker 
lines from singly-ionized iron-group elements are also 
seen. Modeling of the absorption line spectrum shows 
8300<T off <9300K (A2-A41a+, with T cS -8500K pre- 
ferred). This is strongly inconsistent with the observed 
He I emission and Lya and Ly/3 fluorescence lines of O I 
and Fe II, implying the presence of an unseen source of 
ionizing photons. 

Photospheric pulsations are apparent in the 
metal absorption lines, with periodic broadening and 
the development of excess blue- wing absorption sugges- 
tive of periods of increased mass loss. However, no evi- 
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dence is seen for eruptive mass loss, and W243 appears 
to be in a quiescent cool-phase state. 

4.2. The current state of W243 

With L„=7.3xl0 5 L Q (assuming d=4.5kpc; Clark et al. 
2005) and T cff ~8500K, W243 lies near the post-Main 
Sequence track of a ~40M Q star (e.g. Schaller et al. 1992), 
falling as expected between the ~40-55M Q progenitors of 
the Wdl Wolf-Rayet population (Crowther et al. 2006) 
and the <3OM ~07-8V stars at the main sequence turn- 
off (Negueruela et al. 2009, in prep.). W243 appeared on 
the blue side of the HRD as a luminous, early-B supergiant 
when observed by Borgman et al. (1970), Lockwood (1974) 
and Westerlund (1987) before undergoing an (unobserved) 
event that led to W243 moving redwards to the 'yellow 
void' and taking on its current appearance as an early-A 
hypergiant. The final stages of this evolution were seen in 
the early observations described in Paper I, with the star 
remaining generally quiescent in its cool-phase state since 
2003. Our modeling shows Nitrogen to be highly overabun- 
dant, with the depletion of Carbon and Oxygen suggesting 
we are observing significant quantities of CNO-processed 
material. Other a-process elements (Ca, Mg, Si) are also 
enhanced. The implication is therefore that W243 is ei- 
ther in an advanced prc-RSG LBV phase, having already 
ejected its less processed, Hydrogen-rich outer layers, or 
has evolved through the RSG phase and returned to the 
blue side of the HRD, with the current CNO-rich material 
dredged up during the RSG phase. W243 does not display 
the extended circumstellar ejecta expected for either evo- 
lutionary stage, but, as commented on in Paper I, this is 
likely to be rapidly dispersed in the extreme environment 
of Wdl. While the N/O ratio in W243 is similar to AG Car 
(Groh et al. 2009), the latter object is both more mas- 
sive and (considerably) more luminous than W243 and the 
same evolutionary phase cannot be assumed, while post- 
RSG objects such as IRC +10 420 are also overabundant in 
Nitrogen (Oudmaijer et al. 2008) 10 . However, future abun- 
dance studies of the unique coeval population of early-B su- 
pergiants, YHGs and BIa + /WNVL Wolf-Rayet precursors 
found in Wdl offer the opportunity to accurately constrain 
both the evolutionary state of W243 itself and the general 
passage of massive stars as they leave the main sequence 
and evolve through the 'zoo' of transitional objects towards 
the Wolf-Rayet phase. 

4.3. Comparison with other stars 
4.3.1. Cool-phase LBVs 

The visual and near-IR spectrum of W243 shows many 
similarities with other cool-phase LBVs. The Si II dou- 
blet is prominent in cool supcrgiants (Davies et al. 2005), 
and is therefore reported in many objects, while the strong 
near-IR N I lines are also seen at ~B8Ia and later, e.g. in 
the spectra of S Dor (Munari et al. 2009) and HR Car 11 
(Machado et al. 2002). The R-band spectrum of W243 



Sodium abundances are used to infer a post-RSG phase for 
p Cas (e.g. de Jager & Nieuwenhuijzen 1997), but the only fea- 
tures apparent in our spectra result from interstellar absorption. 
11 The N I lines are notably stronger in W243 than in HR Car, 
implying a later spectral type (Munari & Tomasella 1999). 



shows strong similarities with the A2Ia+ LBV HD 160529 
(Stahl et al. 2003; Chentsov et al. 2003), although subtle 
differences are present, with the C II AA6578, 6583 dou- 
blet notable in the spectrum of HD 160529 and Ne I also 
reported; neither can be robustly identified in the case of 
W243. The rich Fe II spectrum is also a typical feature of 
cool-phase LBVs, with HR Car a canonical example, dis- 
playing prominent Fe II emission/absorption blends in the 
cool states of 1991 and 1998-2002 that fade as the LBV 
moves bluewards on the HRD. The transition from P Cygni 
to inverse P Cygni profiles in Fe II seen in Figure 8 is also 
reported in the cases of S Dor (Wolf & Stahl 1990; Wolf 
1992), R127 (Wolf 1992; Walborn et al. 2008) and HR Car 
(Nota et al. 1997; Machado et al. 2002; Crowther, priv. 
comm. 2009); we return to this issue in Section 4.5. 

However, notable differences also exist between W243 
and other early-A LBVs. The Ha and H/3 lines seen in the 
spectrum of W243 are unusual, appearing strongly in emis- 
sion but lacking the P Cygni absorption component gener- 
ally seen in LBV spectra. The strong He I and Lya-pumped 
Fe II and O I emission in W243 are also clearly discrepant 
with the cool-phase emission and absorption lines and are 
more in keeping with early-B/Ofpe hot- phase LBVs (e.g. 
AG Car; Groh et al. 2009), although other characteristic 
hot-phase features such as emission lines of Fe III, N II, 
Si II and Si III are absent. The He I AA5876, 6678 lines 
are seen in absorption in the cool-phase spectra of HR Car 
(Machado et al. 2002), and the four B6-A2 (candidate) 
LBVs described by Chentsov et al. (2003), while He I lines 
appear in emission in R127 only as the Fe II spectrum fades 
during the transition to the hot, quiescent state (Walborn 
et al. 2008). HD 160529 provides a striking example of 
these differences, with the otherwise-similar A2Ia + spec- 
trum showing He I in absorption and a strong P Cygni 
profile in the Ha and H/3 lines; our model also suggests that 
He I should be visible weakly in absorption, rather than the 
observed strong emission. Infra-red spectra of the mid-B or 
later (candidate) LBVs listed by Groh et al. (2007) again 
show W243 to be atypical, with a strong He I A1.083^m 
emission line (of slightly greater strength than the adja- 
cent Pa7 line) contrasting with the remainder of the sam- 
ple that show the line to be in absorption in two cases 
(HD 168625, HD 168607), almost absent in HR Car and 
showing a weak emission component with strong P Cygni 
absorption in HD 160529. 

4.3.2. Yellow Hypergiants 

Our modeling shows that the absorption line spectrum of 
W243 is well described by a cool hypergiant close to the 
Eddington limit, and we observe spectral variability that 
is similar to the well-studied YHGs p Cas (Lobel et al. 
1998) and HR 8752 (dc Jager & Nieuwenhuijzen 1997). 
Within Wdl, W265 also displays pulsational variability in 
the near-IR N I lines, while at radio wavelengths W243 
is very similar to the F2Ia+ YHG W4 (Dougherty et al. 
2009)'. The lack of Fe I features indicates that W243 is hot- 
ter than these objects, which show strong neutral metal 
absorption and a complete absence of He I, with the latter 
also not observed in the A- type YHGs W12a (A5Ia + ) and 
W16a (A2Ia+). The peculiar mid-A YHG IRC +10 420 
(Oudmaijer 1998; Humphreys et al. 2002) shows the 
strongest spectral similarities to W243, with Fe II multi- 
plet 40 and 74 lines in emission and a strong, double-peaked 
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Ha emission line with broad electron-scattering wings; both 
objects also display double-peaked Ca II emission. However, 
the close agreement between the Ha and Ca II features in 
velocity space (e.g. Figure 5, Humphreys ct al. 2002) is 
not seen in W243 12 , suggesting that the origin may not 
be the same. In addition, the 2009 VLT/FLAMES spectra 
show that the Ca II infra-red triplet no longer displays a 
double-peaked profile, with only a weak 'step' in the red 
flank of the emission line hinting at the previous location 
of the absorption feature. Lack of contemporaneous R-band 
spectra means that we cannot tell if the double-peaked Ha 
profile is also now absent, although we note the clear weak- 
ening in the redwards peak in the VLT/UVES spectra and 
its apparent absence in the lower-resolution VLT/FLAMES 
LR6 spectra from 2005 (see Figure 2). Therefore, although 
the double-peaked line profiles seen in W243 may corre- 
spond to the disk (e.g. Jones et al. 1993) or bipolar out- 
flow (Oudmaijer ct al. 1994; Davies et al. 2007) models 
proposed for IRC +10 420 13 , the central absorption may 
also result from radiative transfer effects in the line-forming 
region; additional modeling would be required to confirm 
this. 

4.4. Origin of the emission line spectrum 

While the spectrum of W243 shows similarities with both 
cool, A-type LBVs and A- and F-type YHGs, the juxtaposi- 
tion of a cool hypergiant absorption line spectrum with hot- 
phase LBV He I and Lyman-pumped Fe II and O I emission 
lines is an obvious discrepancy. However, our model is un- 
able to reproduce any of the observed emission features 
while simultaneously preserving the N I and Si II absorp- 
tion lines, and we interpret this as an indication that W243 
is a binary, consisting of the cool-phase LBV and an unde- 
tected hot OB (or possibly Wolf-Rayet) companion that is 
the source of the ionizing flux responsible for the observed 
radio emission and discrepant emission lines. The binary 
fraction amongst massive stars in Wd 1 is known to be 
large (Crowther et al. 2006; Clark ct al. 2008; Ritchie et 
al. 2009), and although W243 lacks the strong non-thermal 
radio and X-ray emission that provide evidence for colliding 
winds in a binary system, weak X-ray emission is observed 
(Clark et al. 2008). Despite the implied presence of a hot 
companion to the LBV, radial velocity measurements are 
not obviously consistent with binarity unless we are view- 
ing W243 at a highly favourable angle or the companion is 
in a long-period, possibly highly-eccentric orbit around the 
LBV primary as has been suggested for rj Carinae (Nielsen 
et al. 2007). An alternative configuration may be similar 
to that observed in the YHG HR 8752, where the flux from 
a B1V companion in a wide orbit is responsible for the 
formation of the observed [N II] emission and a compact 
radio nebula, neither of which could be formed by the cool 
hypergiant itself (Stickland & Harmer 1978). A similar sce- 
nario may also apply to the YHG W265, which also displays 
[N II] and radio emission (Ritchie et al. 2009). In the case 



This may be a result of the blending with the adjacent 
Paschen-series absorption/emission features which are not seen 
strongly in IRC +10 420 

13 Humphreys et al. (2002) also propose a complex rain model 
for IRC +10 420 in which both inflow and outflow are present 
in an opaque wind. However, this requires a mass-loss rate ~10 3 
times greater than currently appears to be the case for W243. 



of W243 a hotter companion is probably required to pro- 
vide the ionizing flux necessary to form the observed He I 
and Lya-pumped emission, both absent in HR 8752 and 
W265, although we note the absence of higher-excitation 
emission lines such as N II, [S III] , [Fe III] and [Ar III] that 
are observed in the sgB[c] binary W9 (Clark et al. 2005, 
2008), which also displays far stronger Lya-pumped lines 
than those observed in W243. 

The lack of observed Fe I or Fe III features in the spec- 
trum of W243 is unusual and implies that almost all of the 
iron is in the Fe + state in the line formation region. This 
is reminiscent of the B and D Weigelt blobs of -q Carinae 
(Verner et al. 2002), which display spectra dominated by 
Fe II and [Fe II] emission originating from upper levels 
around 6eV (permitted lines) and 2eV (forbidden lines). 
This Fe II spectrum arises from collisional excitation at 
10 6 < n e < 10 7 cm~ 3 populating the lower Fe II levels, 
with strong continuum pumping routes to ^6eV levels aris- 
ing from the a 4 D and a 2 G levels (the a 2 G level is the upper 
level of our observed [Fe II] multiplet 14F emission) and 
Lya pumping from the a 4 D and a 4 G levels; this also ap- 
pears a plausible model for the Fe II emission region of 
W243 14 . The significant variability in Fe II, He I and O I 
emission suggests that the efficiency of the pumping chan- 
nels varies, possibly due to changes in Lyman-series optical 
depth or other radiative transfer effects. 

4.5. Location of the emission line formation region 

The transition from P Cygni to inverse P Cygni profiles in 
the Fe II lines seen in Figure 8 is observed in other LBVs, 
and is interpreted as an indication of infall during the ra- 
dial contraction phase of the stellar photosphere during the 
S Dor variability cycle (Wolf & Stahl 1990). In the case of 
W243, this transition occurs at a time when the strong Si II 
and N I absorption lines develop broad blue wings 15 but in- 
spection of the Fe II multiplet 74 and higher Paschen-series 
lines suggests that these may not be true (inverse) P Cygni 
profiles, and rather an effect of emission at near-constant 
radial velocity superimposed on an absorption line with a 
varying line centre: this is similar, but more pronounced, to 
the effects seen in p Cas where static emission lines and vari- 
able absorption lines produce apparent line-splitting (see 
Lobel et al. 1998, figs. 13 and 15). With W243, if the ab- 
sorption is at its most blucshiftcd (and especially if a broad 
blue wing also develops) the emission/ absorption blend will 
take on an apparent P Cygni profile (e.g. as visible in the 
Pa<5 line in the left panel of Figure 2 in the 10/07/2004, 
MJD=53192.2 spectrum, or in the Fe II multiplet 40 lines 
in Figure 8), while at its most redshifted the absorption 
component will be seen entirely on the red flank of the 
emission component, creating an inverse P Cygni profile: 
however, neither case corresponds to the significant outflow 
or inflow normally associated with these spectral features. 



1 The similarity between the Fe II-rich spectra of W243 and 
other cool-phase LBVs, e.g. HR Car or S Dor, also suggests 
similar conditions in the Fe II line formation regions, although 
we do not suggest that a hot companion is required in all of these 
objects; instead, the similarity probably reflects the complex 
balance of collisional, continuum and Lya-pumping processes in 
these objects (Verner et al. 2002). 

15 A similar transformation is correlated with periods of en- 
hanced mass- loss in the YHG p Cas (Lobel et al. 2003) . 
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The near-static radial velocities of the undistorted emis- 
sion lines over the course of our observations and the close 
correspondence between the radial velocities of the permit- 
ted and forbidden lines suggests that they all form at sig- 
nificant radii, when material expanding slowly at constant 
velocity is ionized by the flux from the hot secondary. As an 
emission spectrum was not observed by Westerlund (1987), 
it appears that this material is a relic of mass-loss during 
the transition from an early-B supergiant to the current 
cool state. The O I A8446, Fe II and higher Paschen-series 
lines have very similar profiles, and all likely form in a com- 
mon region in the boundary zone between H I and H II 
regions required for efficient Ly/3-pumping of the O I line, 
while the consistency between the measured radial veloci- 
ties of almost all the permitted and forbidden emission lines 
suggests formation in the same comoving region. The radial 
velocity of the Ca II emission is somewhat variable, while 
the line is broader than all but the Balmer-series emission 
lines, most probably reflecting formation in a different lo- 
cation: the ~11.9eV ionization potential means recombina- 
tion will be negligible in a H II region, and the line likely 
results from resonance absorption in an H I region rather 
than being limited to the H I/H II transition zone. It is 
plausible that the Ca II emission may be a result of dust 
evaporation (e.g. Prieto et al. 2008) which would indicate 
a post-RSG state for W243, although the [Ca II] emission 
that would provide stronger support for this hypothesis is 
absent (Shields et al. 1999). These lines lie in a region of 
heavy telluric contamination, but we also note that [Ca II] 
is only weakly detected in the B[e] star MWC 349 A despite 
strong Ca II emission being present (Andrillat et al. 1996). 
The inferred electron density in the line formation region 
is high (Hamann & Simon 1988), and it is likely that the 
Ca II 3 2 D level is collisionally de-excited before the forbid- 
den transition can occur. This also appears likely in the case 
of W243, where relatively strong [Ca II] emission would be 
required to be detectable amidst the telluric lines. 

4. 6. Conclusions and future work 

Our observations of W243 show that it has remained in a 
quiescent, cool-phase state since 2003, with no indication 
of either further movement redwards or a return to its for- 
mer hot state. The ~A3Ia + absorption-line spectrum dis- 
plays photosphcric pulsations and episodic mass loss that 
appears similar to warm YHGs (e.g. p Cas, Lobel et al. 
1998 or W265, Ritchie et al. 2009). The high Nitrogen 
abundance and depletion of Oxygen and Carbon imply that 
CNO-processed material is on the surface of W243, indicat- 
ing an advanced evolutionary state and significant previ- 
ous mass-loss. Superimposed on the YHG absorption-line 
spectrum are emission lines of H, He I, Ca II, Fe II and 
O I formed at large radii when material lost from the LBV 
is ionized by an unseen hot companion. We are unable to 
strongly constrain the emission-line formation region and 
the origin of the peculiar, double-peaked Balmer-series pro- 
files or the highly variable He I emission lines. However, 
high-resolution infra-red spectroscopy with adaptive optics 
offers the opportunity to directly probe the line-forming re- 
gion and thereby determine the recent mass-loss history of 
W243. In addition, abundance studies of the coeval popu- 
lation of both pre- and post-LBV transitional stars in Wdl 
will allow the evolutionary state of both W243 and its evo- 
lutionary contemporaries to be determined. 
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Appendix A: List of identified absorption lines in 
the spectrum of W243 

W243 displays a large number of absorption lines from neu- 
tral and singly-ionized metals; these are listed by species 
and multiplct in Table A.l. To identify the lines wc fol- 
lowed the approach of Oudmaijer (1998), first identifying 
strong absorption lines and then searching for other lines 
in the same multiplct using the line identification tables of 
Moore (1945). We then searched for other multiplets of 
the same species, and regard as robust any identification 
where several lines from the same multiplet are unambigu- 
ously detected. In a few cases tentative identifications are 
made when only one line in a given multiplet is detected, 
provided that one of the following conditions apply: 

— the multiplet contains no other lines, or the other ex- 
pected lines in the multiplct arc outside the range of 
coverage of the spectrum. 

— the proposed identification is from a species that has 
been identified in other multiplets, the line is not close 
to possible interstellar or telluric features and no other 
plausible identification exists. 

— the proposed identification is from an otherwise uniden- 
tified species, but the line is close to a robustly identi- 
fied line (allowing the difference between line centres 
(5 A = Aknown — A un known to be measured accurately) and 
no other plausible identification exists with the mea- 
sured SX. Nevertheless, such an identification should be 
considered tentative. 
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The majority of absorption lines in the range 4800 
to ~7000A could be positively identified in this manner. 
After eliminating obvious telluric features, many remain- 
ing absorption features could be identified as being due to 
diffuse interstellar bands by comparison with the Diffuse 
Interstellar Band Catalog 16 , which is based on Jenniskcns 
& Desert (1994). The large number of telluric features 
above 7000A made comprehensive identification in that re- 
gion problematic, and our list is incomplete in that wave- 
length range with only the stronger features being identi- 
fied. 



http:/ /leonid. arc. nasa.gov/DIBcatalog. html 
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A.l. Identified absorption lines in the spectrum of W243 between 4800-9000A. 
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Table A.l. continued. 
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^Multiplet from Moore (1945). Lines that are unclassified are marked as 'unci.'. 
a inverse P Cygni profile in 21/09/2003 spectrum 
fc blended with nearby line 
c uncertain 

^highly uncertain, very weak 

e excess absorption apparant on blue wing 

^ possible P Cygni profile 



